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Abstract. Although known for almost a century, the photophoretic force has only recently been considered in astrophysical 
context for the first time. In our work, we have examined the effect of photophoresis, acting together with stellar gravity, 
radiation pressure, and gas drag, on the evolution of solids in transitional circumstellar disks. We have applied our calculations to 
four different systems: the disks of HR 4796A and HD 141569A, which are several Myr-old AB-type stars, and two hypothetical 
systems that correspond to the solar nebula after disk dispersal has progressed sufficiently for the disk to become optically thin. 
Our results suggest that solid objects migrate inward or outward, until they reach a certain size-dependent stability distance 
from the star. The larger the bodies, the closer to the star they tend to accumulate. Photophoresis increases the stability radii, 
moving objects to larger distances. What is more, photophoresis may cause formation of a belt of objects, but only in a certain 
range of sizes and only around low-luminosity stars. The effects of photophoresis are noticeable in the size range from several 
micrometers to several centimeters (for older transitional disks) or even several meters (for younger, more gaseous, ones). We 
argue that due to gas damping, rotation does not substantially inhibit photophoresis. 
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1. Introduction 



The standard (core accretion) scenario of planet formation 
implies continuous growth of soli ds in a protoplanetary disk 
aroun d a young star (s e e, e.g. ISafronovl 1 1 972k IWetherilj 



1993 



198(1 IShuetalJ 119871 iLissauetj 11993b IStro m & Edwards 



Weidenschilling & Cuzzilll993t 
Blum & Wurml l200rj IWurm & Blum 



Papaloizou et al. 1999; 



2006: iHenningetal 



2006; M eyer et alJ 120071) . Micron-sized dust grows step by 
step to larger bodies, until asteroid-sized planetesimals are 
reached. At sizes well below 1 km, gravitational interac- 
tions between the objects play a minor role and the pro- 
cess is largely determined by interactions of solids with the 
ambient gas in the disk. Gas causes sedimentation, mix- 



ins, radial drift dWeidenschilling 


1977; 


Klevetal. 1993: 


Takeuchi & Artymowicz 2001; Brauer et al. 


2007) and other 



effects that all dictate the spatial, size, and velocity distribu- 
tions of grains and therefore the conditions for their growth. 
As the disk evolves, both gas and dust are gradually removed 
from the systems. In this process, the gas-to-dust ratio re- 
duces from ~ 100 in the initial protoplanetary disk stage, at 
the ages on the order of 1 Myr, to vanis hing values of <C 1 
at the debris disk s t age after ~ 10 Myr dLawson et al.ll2004l 
Haisch et al] 120051; iHollenbach et all 120051; iTakeuchi et alJ 
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2005 1 : IJavawardhana et al]2006l:lBalog et alj2007l;ICurrie et al 



20071) . 

Transition from gas- and dust-rich, optic ally- thick pro- 
toplanetary disks to nearly gas-free, optically-thin debris 
disks is currently in the f ocus of interest of both observa- 

1 1 «-i — l i\\K T + — Hi r\o?l li/r~-f:iil 



1983 




1988; Weidenschilling & Cuzzi 


1993; Strom et al. 


1993 


Zuckerman & Becklin 1993; Simon & Pratol 


1995; 


Mever & Beckwith 2000; Ardila et alj2005l;ICalvet et alj 


2005 


Hueso & Guillot 2005; Augereau 2006; Bouwman et al. 


2006; 


Eisner et al.1 120061). Examples of transitional objects are TW 



Hya, HR 4796A and HD 141569A, while the well-known sys- 
tem j3 Pic already enters into the realm of almost gas-free de- 
bris disks dThebault & Augereaull2005l) . Resolved transitional 
disks exhibit radial structure in the form of rings and gaps, of- 
ten alternating, and there is a lot of debate whe ther this struc- 
ture is caused by gravity of hidden planets dAugereau et al] 



1999b; IWvattet all 1 19991; iTelesco et al.l BOQOb or by interac 



tion of solids with the ambient gas component dKlahr &Lin 
2000; ITakeuchi & ArtvmowiczfeoOll: iBesla & Wull2007t) . The 



first possibility is supported by the fact that the amount of gas 
remaining in these systems i s probably no longer sufficient to 
form Jupiter's gas envelope dChen & Kampll2004l) and so, the 
planet formation process must already be finished. The sec- 
ond hypothesis is substantiated by simulations based on the 
observational estimates of the gas contents. These show that 
gas drag, acting together with other forces - stellar gravity and 
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radiation pressure — could result in segregation of different- 
sized solids and radial fractionation of dust. 

In this paper, we study the motion of different-sized solids 
under the combined action of stellar gravity, radiation pressure, 
and gas drag, to which we add photophoresis, an additional 
radial force that acts on particles in a gas disk which are ex- 
posed to the radiation field of the central star. Photophoresis, 
one of the lesser known for ces of physics, has been first de- 
scribed by lEhrenhaftl dl917l) . The force is caused by the fol- 
lowing process. If an object is embedded in a thin gas and 
is exposed to an anisotropic radiation field (i.e. of a star), a 
gas molecule that becomes accommodated by the object's sur- 
face and rejected again will depart from the illuminated - and 
therefore warmer - side on average with a greater velocity 
then from the dark, colder one. Net momentum is transferred 
to the object, accelerating it away from the light source (note 
that for very small partic les, the force can also be attractive - 
Tehranian et al.ll200lb . Although known for almost a hun- 



see 

dred years and successfully used in tec hnical applications, such 
as the construction of optical traps dSteinbach et al.1 120041) . it 
has only recently been analyzed in astrophysical context. The 
first attempt t o examine its influence on planetary formation 
was made bylR rauss & WurrrJ (120051) . In their subsequent pa- 
per (^^rm|&&aus3200g), they investigated how photophore- 
sis might affect the formation of chondrules and the survival 
of Ca lcium/ Aluminum-rich inclusions. Finally, iRrauss et al.1 
(2006) analyzed the effects imposed by photophoresis on the 
inner rim of a dusty protoplanetary disk. They assumed the disk 
to stay optically thick at all times, thus allowing only its in- 
ner edge to be influenced by photophoresis. Using an exponen- 
tial la w to describe particle grow th and t he a 1+1D tu r bulen t 
model (Shakura & Sunvaev 1973) used in Alibert et al. (2005) 
and lPapaloizou & Terqueml d 1999b for the evolution of the gas 
disk, Krauss and Wurm computed the motion of the inner edge. 
This is identical to the motion of the smallest particles, since, 
as will be shown subsequently, for solar luminosity systems, 
the velocity of outward motion of an object increases with its 
size in a certain size range (compare Fig. |3). In essence, they 
showed that the dust disk's edge first moves outward, pushed 
by photophoresis. After a few million years of time, it shrinks 
towards the star again, because gas pressure declines as disk 
dispersal proceeds. 

This paper focuses on transitional disks. In these systems, 
the conditions are quite favorable for photophoresis to be effi- 
cient. Transitional disks are already optically thin, so that the 
particles are exposed to the strong radiation field. On the other 
hand, they are still sufficiently gas-rich. We would like to check 
whether, and to which extent, photophoresis may affect the ra- 
dial migration of solids. 

We use a single-body dynamics approach to investigate 
the behavior of a two-dimensional swarm of particles, the 
sizes of which are distributed in an interval ranging from 
several micrometers to several meters. The size of an indi- 
vidual particle, though, is being kept constant. We develop 
a simple theory for the long-term behavior of a distribu- 
tion of solids experiencing gravity, radiation pressure, gas 
drag and photophoresis. Metho dically, our study is similar to 
Takeuchi & Artvmowiczl d200ll) . whose computations we gen- 



eralize to the presence of photophoresis. To solve the equa- 
tion of motion, we employ a mod i ficatio n of the semianalytic 
scheme used by Weidenschillinj ( 1977 ). The results are ap- 
plied to circumstellar disks with different gas density around 
stars of different luminosity. 

In Sect. 2, we discuss the astrophysical status and parame- 
ters of the systems chosen for this study, as well as the formu- 
las used for gas density and temperature in the disks. Section 3 
provides formulas for different perturbing forces, including the 
photophoretic one. Section 4 deals with the equation of motion. 
In Sect. 5, we analyze the properties of the particles' radial mo- 
tion and derive a formula for the radii of stable orbits. We also 
estimate the size ranges in which photophoresis should be taken 
into account. In Sect. 6, we examine the problem of particle ro- 
tation in order to check whether it might nullify the effect of 
photophoresis. In Sect. 7, the results are summarized and ideas 
for future research are presented. 

2. The systems 

Following iTakeuchi & Artvmowicz (2001), we choose two 



transitional disks of particular interest for our study: HR 4796 A 
and HD 141569A. We use these two systems in our work as 
model environments for the forces we wish to explore, as their 
properties are quite well known and reasonably accurate mod- 
els exist for the physical conditions (gas density and tempera- 
ture) in their disks. 

The AO-type star HR 4796A, located at a distance 
of 67.II34PC from the sun (Hipparcos data), is a mem- 
ber of t he TW-Hydrae-As sociation (TWA in the follow- 



ing, see iKastner et alj 1200 lb . According to current knowl- 
edge, this is the young s tar association closest to our sun 
dZuck erman & Songll2004 ). The age of TWA can be estimated 
to be approximately 10 Myr. Since it has used up its molecular 
cloud completely, the only gas left in in the association is bound 



L + = 21.0L m ( 


Koerner et al. 


Telesco et al. 


2000). IRAS 



are es t imated to be = 2 . 5M ( 
11998k IJayawardhana et all \l99 
discovered that it emits 0.5 % of its entire radiative power in the 
infrared part of the spectrum, which points to a high dust den- 
sity in its neighborhood - in fac t, HR 4796A is the du st-richest 
star in the Bright Star Catalog. St auffer et al.1 d 1995b estimate 
its age to be 8 ± 2 Myr, which makes it slightly younger than j3 
Pic (12 Myr). The structure of the dust disk around HR 4796 A 
is highly complicat ed. First resolved im age s in thermal infrared 
were o b tained bv | Koerner et al. ( 1998b and lJavawardhana et al 



dl998b . IWahhai et al.l d2005b conducted elaborate studies of 



the disk's structure, using data from MIRLIN at the Keck II 
telescope as well as 350 jum measurements from the Caltech 
Submillimeter Observatory and Hubble Space Telescope scat- 
tered light images. Their studies point to a disk composed of an 
inner, exozodiacal dust ring located at roughly r = 4AU from 
the star and a wide, two-component outer dust belt consisting 
of a broad ring of ~ 7^im grains stretching from 45 to 125 
AU, and a narrower structure between 66 and 80 AU consist- 
ing of ~ 50/im grains. While the exozodiacal dust may be the 
product of an asteroid-type belt, the outer belt can be explained 
most naturally by the assumption that the grains are emitted by 
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Table 1. Parameters of gas disk models. 



System 


Example 


p„[10- lu kgm- J ] 


Density exponent q 


Luminosity L*/L© 


Mass M^/Mq 


1G1L 


'older' SN 


3.12 


-2.75 


1.0 


1.0 


hGIL 


'younger' SN 


156 


-2.75 


1.0 


1.0 


lGhL 


HR 4796A* 


3.12 


-2.75 


21.0 


2.5 


hGhL 


HD 141569A* 


156 


-2.25 


22.4 


2.3 



*Sources: Takeuchi & ArtvmowicJ d200lh . their Table 1 



an exo-Kuiper swarm of planetesimals. Since the smallest ones 
are diffused rapidly by radiation pressure, they form a broad 
belt, while the larger ones stay closer to their a rea of origin 



For gas density and temperature in all four model systems, 



Also, asymmetric structures have been reported dTelesco et al 



2000). They may be caused by secular perturbations of one or 
several planets or by the stellar companion HR 4796B. Most 
probably, the disk is perturbed by planetary bodies as well as 
by the B component. In order to define the mass and number of 
planets, the orbit of HR 4796B would have to be known with 
greater precision. 



The second object, HD 141569A, has 
distance of 99 ± lOpc. Its par ameters 
as M+ = 2.3M W , L* = 22AL R djura et al 



the Hipparcos 
are estimated 



1993, 1998 



van den Ancker et al . 1998) It is underluminous for its spectral 
type B9.5 Ve, which is a common occurrence for you ng A-type 
stars, found also at HR 479 6A, j3 Pic and 49 Ceti djura et al 



lLowrance etal]|2000h . Being, like HR 4796A, part of 
a multiple system with an M2 and an M4 component, the 
lower-mass companions which have not yet reached the main 
sequence c an be used to find an esti mate for the system's age of 
5 ± 3 Myr dWeinberger et al.ll2000h . As in th e HR 4796A case, 
the d i sk itself has a complex morphology (We inberger et al 
1999 1: lAugereau et alJll999aL iFisher et al]l2000tlMouillet et al 



200 lb . consisting of two dust belts with radii of 200 and 325 
AU, the centers of which are shifted by 20 - 30 AU in the 
direction of the system's semiminor axis. Between them a 
dust-free gap can be found at ~ 250 AU. Also, inside of 150 
AU the disk's luminosity declines to the level of background 
noise, which points to a strong depletion of dust in the system's 
inner region. The outer ring shows a tightly-wrapped spiral 
struc ture, which, according to t he numerical calculations 
of lAugereau & Papaloizoul d2004l) . can be produced by the 
gravitative perturbation of HD 141569B and C. The dust gap 
at 250 AU may be caused by a planet of approximately Jovian 
mass, but it is not yet clear whether gas giants can form within 
a few mi llion years a t a distance of several hundred AU from 
the star (Wyatt 2005|). 



Both disks - HR 4796 A and HD 14 1569 A - surround lumi- 
nous, early-type stars. Since we wish to explore the dependence 
of the particle dynamics on the stellar luminosity as well, we 
add two more, hypothetical systems with gas densities equal 
to those of HR 4796A and HD 141569 A, but around a star of 
solar luminosity. They can be viewed as repr esentations of the 



solar nebula at different stage s of dissipation (Hol lenbach et al 
1994tlHollenbach et alJEoOOl) . 



we u se standard power-law approximations dHayashi et al 
1985): 



1/4 



(lAu) 



-1/2 



and 



lkgm 



-3 



PQ 



1 AU 



(1) 



(2) 



where po and q are constants that vary from one system to an- 
other. To keep the treatment simple, we assume the exponent q 
to be constant throughout each disk and do not use a 'break-off 
func tion' to describe the density dr op beyond a certain radius 
r out dTakeuchi & Artvmowical200ll) . As will be shown subse- 
quently, the outer part of the radial profile does not really mat- 
ter, because photophoresis affects only the inner parts of the 
disk. 

Altogether, our model contains three basic parameters: the 
gas density po at r = 1 AU, the gas density exponent q which 
controls the size of the gas disk, and the stellar luminosity L*. 
Our four model systems essentially explore the photophoretic 
effect in the 'parameter rectangle' luminosity - gas density: 

- ILhG (low luminosity, high gas content, corresponding to 
the solar nebula at an earlier stage), 

- ILIG (low luminosity, low gas content, solar nebula at a later 

stage), 

- hLhG (high luminosity, high gas content, HD 141569A), 
and 

- hLlG (high luminosity, low gas content, HR 4796A). 
The parameters of all four systems are listed in Table Q] 

3. Forces 

In circumstellar nebulae, solid bodies experience a number of 
perturbing forces in addition to gravity, causing them to move 
along non-Keplerian orbits. In a first approximation, we can as- 
sume that they describe circular orbits, the radii of which shrink 
or grow in the course of time, depending on the size of the body 
and its distance from the star. We now introduce the different 
forces, presenting an analytic expression for each. 

3.1. Photophoresis 

While Rohatschekl d 1996b found a semi-empiri cal formula for 



the ph otophoretic force at all gas pressures, iBeresnev et al 



( 1993) derived an analytic expression for spherical, nonrotating 
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objects with a homogenous surface using an elaborate theoret- 
ical approach starting f rom the molecular v elocity distribution 
function (compare also lTehranian et al .120011) . The acceleration 
due to the photophoretic force they found is: 



2phot 



4*Pbuik 



2kT 



a E + 15AKn (1 - a E )/4 + a E A x i> 2 



(3) 



Here, s is the particle radius, / the radiation intensity, Ji the 
asymmetry parameter that describes the accommodation of 
gas molecules to the particle's surface and light absorption. 
Assuming complete absorption and an accommodation prob- 
ability of 100 %, we set J\ = 0.5. The mean molecular weight 
of a gas of solar composition is denoted by — 2.34 and the 
mass of a hydrogen atom by ran- Further, k is Boltzmann's con- 
stant, T the gas temperature and Pbuik the density of the solid 
material. For icy aggregates, Pbuik ~ 1000 kgm~ 3 is a good ap- 
proximation. The energy accommodation coefficient a E is the 
fraction of molecules in contact with the surface that accommo- 
date to the local temperature, which enables them to contribute 
to photophoresis. In our work we assume complete accommo- 
dation, and thus set a E = 1 . The heat exchange parameter A, 
which describes the particle's thermal relaxation properties, is 
determined as 



A = 



A: th + 4ec7r 3 5 



(4) 



tgas 



where is the material's thermal conductivity, e its emissiv- 
ity, a the Stefan-Boltzmann constant and £ gas the thermal con- 
ductivity of the gas. Assuming black bodies, we set e = 1. The 
quantities *Pi and ^2 are given by 



¥1 



Kn 



2% x l 2 Kn 



Kn + (57T/18) 



' 1 2 4 



5Kn 2 + % l l 2 Kn + n/4 
\.2\n x > 2 Kn 




lOOKn 2 + ic/4 



(5) 



and Kn is the Knudsen number. It is defined as Kn = L/s, where 
L = l/( \/(32) 7rn gas r 2 ) is the molecule's mean free path, with 
r g w 10~ 10 m being the radius of gas molecules and n gas their 
number density. 

Expression Eq. (01 for the photophoretic acceleration is 
valid for all Knudsen numbers. In the case of high Knudsen 
numbers (mean free paths of molecul es are large compared t o 
particle sizes), Eq. (01 reduces t o (seelKrauss & Wurml d2005h . 
their Eq. 1, and lBeresnev et all ( 1993 ). their Eq. 26) 



-*phot 



where 



4pbulk (The + Tad + T gas ) 



The = kfoT , 

T rad =4aT 4 es, 
T gas = p J 2kT/7tH g mu s . 



(6) 



(7) 
(8) 

(9) 



hLIG (HR 4796A) 



1e+12 




100 



rl AU 
hLhG (HD 141569A) 




100 



r/AU 



Fig. 1. Knudsen numbers Kn = L/s for systems with low (top) 
and high (bottom) gas content as a function of distance from 
the star. In each panel, four lines correspond to the particle radii 
of 1 mm, 1 cm, 10cm, and 1 m. Since L only depends on gas 
density, Kn does not vary with stellar luminosity. 



Here, p = n g!ls kT is the gas pressure (assuming ideal gas). In 
transitional disks, we can usually assume to be in the high 
Knudsen number regime, except for large objects in the inner- 
most parts of the disks (Fig. [TJi. The simplified expression © 
will allow us to find useful approximate solutions for the parti- 
cles' stability radii (i.e. the radii of stable circular orbits). For 
large particles (s > 10 cm) and small distances from the star 
(r < 1 AU), the more general expression (01 will be used in nu- 
merical calculations. 

The terms Yh c , T ra d and Y gas are related to three different 
processes which reduce the photophoretic effect. The first term, 
Yhc, corresponds to the transport of thermal energy through the 
object - i.e. heat conductivity, which reduces the temperature 
gradient, thus lowering the efficiency of photophoresis. The 
second one, Y ra d, describes thermal radiation from the parti- 
cles' surface. Since it is proportional to T 4 , much more energy 
is radiated from the warm than from the dark side, which re- 
sults in a substantial reduction of the temperature gradient, and, 
therefore, of photophoresis. Obviously, this process gains effi- 
ciency with increasing temperature. Finally, the third term Y gas 
describes heat conduction away from the object's surface into 
the surrounding gas. Since it also shows a weak dependence 
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on temperature (<* T 1 / 2 ), it causes a reduction of temperature 

gra dient too. 

Krauss & Wurrril d2005l) and IWurm & Kraussl d2006h use 



only the first term Th c in the denominator of Eq. (O, which is a 
good approximation for particles with s <C 1 mm and only in the 
case of low luminosity (i.e. solar-type) systems (see Sect. 15.2b . 
In order to analyze the relative importance of the three 'coun- 
tereffects' Yh c , Y ra( j and T gas , we plot them for HD 141569A as 
functions of particle size (Fig. |2j at three different distances. 



1000 
100 
10 

1 

0.1 
0.01 
0.001 
1e-04 
1e-05 
1e-06 




1e-07 



100 



1000 



10000 
s / (im 



Y 

gas 

100000 



1e+06 



Fig. 2. Quantitative comparison of the three reduction pro- 
cesses for hLhG (HD 141569A) at r = 5, 10,20AU (from 
above for each set of lines). For the other systems the result 
looks qualitatively similar, with T gas being in the case of the 
gas-poorer systems completely negligible for all sizes s. 

Figure [2] shows that for fairly small particles (s < 1 mm), 
heat conduction (Th c ) is the main reduction process. At larger 
sizes (s > 10 cm), radiation (Y ra d) takes over. In the intermedi- 
ate interval, both processes contribute to the reduction of pho- 
tophoresis. At yet larger sizes (s > 1 m), heat conduction from 
the surface into the surrounding gas (T gas ) comes into play, ap- 
proaching values similar to that of Th c - As lower gas pressure 
dramatically reduces the importance of T gas , it is completely 
negligible for all object radii in the case of the gas-poor sys- 
tems hLlG (HR 4796A) and 1L1G. 



3.2. Radiation pressure 

Since radiation pressure is proportional to r~ 2 , it can be 
taken into accoun t by introducing an 'effective stellar mass' 
dBurns et alJll979t) : 



M e ff = M*(l-j3), 



(10) 



where j5 stands for the ratio of accelerations due to radiation 
pressure and gravity: 



J3 = 



0.5738 g pr L*/L G 



■0 



^rad 

flgrav Pbuik[gcm- 3 ]i[^m] M*/M 

, 1 

s 



1 /J.m 



(ID 



with 2pr being the radiation pressure efficiency. We will make 
use of the constant B later. The effective 'photogravitational' 
acceleration is then 



■■grav 



eff = —G 



M eff 



(12) 



In our work, we simply set Q pr = 1, assuming the solids to 
be black bodies. The reason is that generalization to non-black 
surfaces requires a change not only in the emissivity, but also in 
the asymmetry parameter J\ . The latter is not known for real- 
istic materials, and its determina tion requires co mplicated nu- 



merical procedures (see, e.g. iMackowski Il989l) . Analytical 



solutions are possible, but only in the cases of hig h or low ab- 
sorption coefficients (see Arnold & Lewittes1ll982 ). 



3.3. Gas drag force 

Unlike the other three forces (gravity, radiation pressure, pho- 
tophoresis), gas drag is not a radial force. Since it results from 
the momentum transferred to the body by molecules impinging 
on it as it moves through the gas, the force vector is antiparallel 
to the relative velocity Av of the dust grain with respect to the 
gas. 

To calculate gas drag, we first need a model for the motion 
of the gas component of the circumstellar nebula. Since the 
pressure gradient supports gas against stellar gravity, it travels 
on circular orbits with a sub-Keplerian speed 



Vg = V K VT -tj, 
and angular velocity 



(13) 



a, = n K \/T^, (14) 

where vk = y/ GM+jr is the Keplerian circular velocity, Hk = 
vk/V = y GM-,, I r 3 the corresponding angular velocity, and 77 
is the ratio of pressure gradient force to gravity: 



'7 



1 dp 



(15) 



The rj ratio can be rewritten as: 



1/ 1.1 / 10 ' ( i-q 



2.34J \L S , 



-1 



1AU 



1/2 



= E 



1AU 



1/2 



(16) 



If the motion of the particles is subsonic (Av <C vt ~ cs 
where cs is the speed of sound - this is always the case in the 
transitional disks), the gas dra g acceleration is given by (see 
Takeuchi & Artvmowic z 200lb: 



a D 



where 



V'x = 



3pg 



vt Av, 



8kT 



1 

1 4 

= ^ X ( v therm) 



3 \7tH g mft / 

is 4/3 times the mean thermal velocity. 



(17) 



(18) 
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The reaction of particles to the gas drag force critically de- 
pends on their size. Small objects adjust their velocity instan- 
taneously, they are swept along with the gas component. Large 
ones react only sluggishly to gas drag, taking a longer time to 
change their velocity substantially. In order to create a quan- 
titative measure for the tendency of solids to be influenced by 
gas drag, the stopping time is introduced. We denote the time 
needed for a particle injected into the gas to be slowed down to 

times its initial velocity by t s = Av / ad- The dimensionless 
stopping parameter T s is defined by: 



4pbulkSV'K 

3p g rv T 



(19) 



The right-hand side approximation holds in the case of sub- 
sonic motion. It renders the stopping time independent of the 
particle's momentary velocity. 

4. Equation of motion 

The equation of motion of a dust particle is: 



d 2 r _ 

— -j — ftgrav, eff + a phot + a D • 



(20) 



Being interested in the radial motion of the bodies, we now 
consider the radial component of the equation of motion d20l . 
In the reference frame corotating with the gas, it takes the form 



d 2 r 
dt 2 



£2 2 T — flgrav.eff + a phot ■ 



(21) 



The acceleration d 2 r/dt 2 = Ag of the particle in the corotating 
frame can be interpreted as the 'residual gravity' . It computes 
to 

Ag= -a g rav,eff + aphot + ^g>-= (j3-r-^-Tj)0|r, (22) 

where % = a phot/<2grav is the photophoresis-to-gravity ratio. 
Particles experience the inward-directed residual acceleration 
Tjil^r and the outward-directed accelerations due to pho- 
tophoresis (x^l^r) and radiation pressure (/3i2^r). 

Assuming now that the solid particle moves at a Keplerian 
circular speed, we can write 



Ag. 



The relative velocity is then approximately 
Ag 



Av = vk — v B ~ — 



2fln 



(23) 



(24) 



We now use Ag to derive an approximation for the ra- 
dial drift velocity of the particles. We can discern two limiting 

cases: 

- Small particles are swept along with the gas. They stay on 
circular orbits, while moving with sub-Keplerian velocity, 
which causes them to experience the residual gravity pull 
Ag. 

- Large particles are decoupled from gas. They move on 
Keplerian orbits, experiencing a head wind that gradually 
reduces their angular momentum. 



In th e first case, the radial d rift velocity v r = dr/df computes to 
(see IWeidenschillingll 1 9771 sect. 4.1, 4.2): 

Vr, small - h Ag = T s (fi + % ~ 1 ) M K . (25) 

In the second case, the orbit decays at a rate 

Ag _ (0+*-i7)rfl K 



r 2 r 

Vr. large = Al' = 

t s v K f s a 



grav 



These two formulas can be combined 
Takeuchi & Artvmowiczll200ll Sect. 3.3): 



(26) 
into (see 

, -vk. (27) 

In Fig. [3] we plot the absolute value of v r over s for r = 
10 AU, both with and without photophoresis. The points where 
v r = are shifted to larger sizes by photophoresis. Note that for 
larger sizes, radial velocities depend almost exclusively on gas 
density since for them gas drag (rj in Eq.|27]l becomes the main 
driving force of radial migration. 
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Fig. 3. The absolute values of the radial velocities at r — 10 AU. 
The minima correspond to the sizes for which v r = 0. 

In what follows, differential equation (f2Tb for r(t) will be 
studied analytically and solved numerically. In the realm of 
large radial velocities (or equivalently, small particles), r(t) 
should be computed from the 'exact' equation of motion d20l ). 
However, our numerical tests have shown Eq. ( f2Tb to be accu- 
rate enough for all particle sizes larger than ^100 jJ.m. 
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5. Results 



5.2. Equilibrium distance 



5.1. Radial motion 

We used Eq. ( f27b to calculate the radial drift of different-sized 
particles. What happens, in short, is that the particles migrate 
inward or outward, until they reach stable circular orbits on 
which gravity, centripetal force (in the particle's own inertial 
system), photophoresis and radiation pressure balance each 
other, thus permitting circular motion (v r = 0), while at the 
same time Av = 0, i.e. the particle travels at the same speed 
as the gas. 
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Fig. 4. Distance as a function of time for particles of four dif- 
ferent radii in the hLlG system (HR 4796A). Dashed lines: 
without photophoresis, solid: with photophoresis. The lower 
border of the plot corresponds to the evaporation limit r eva p = 
0.16AU (defined by black-body equilibrium temperature T = 
1500K). 



As an example, Fig. [4] shows the r(t )-curves for particles 
of four different sizes with and without photophoresis in the 
system hLlG (HR 4796A). The particles start on circular or- 
bits with ro = 5 AU. In fact, initial circularity is not important, 
because elliptic orbits are quickly circularized by the gas drag 
force. The initial value ro does not matter either, because the 
particles migrate towards their equilibrium orbits within several 
thousand years. While small, outward migrating bodies travel 
slowly towards their equilibrium distance r sta b(s), approach- 
ing it asymptotically, the larger, inward migrating ones almost 
'drop' onto their stability orbits, being stopped almost instanta- 
neously. Note, however, that 'instantaneous' here refers to a de- 
celeration process taking several centuries. Therefore, the large 
particles travel essentially on Keplerian orbits, corresponding 
to the second case described in sect. 4. With increasing particle 
size, though, this process takes longer, as larger objects experi- 
ence weaker drag acceleration than smaller ones. Another con- 
clusion from Fig.|4]is that the stability radii are pushed outward 
by photophoresis. 



We now calculate the equilibrium distance which, according to 
Eq. ( T27l >. must satisfy 



(28) 



From this, an implicit expression for the radius r sta b(s) of the 
stable orbit can be derived: 



r s tab(s) 



*grav 



eff( r stab, S) - a p hot(r s tab, s) 



(29) 
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Fig. 5. Stability radii as functions of particle radius s for two 
high-luminosity systems (top) and two low-luminosity ones 
(bottom). In the bottom panel, the 'No Photophoresis' curve 
is identical for both systems, as they differ only in the value 
of gas density, which does not affect the stability condition if 
photophoresis is not taken into account (r\ only depends on the 
density exponent q, not its absolute value - see Eq.[ToT>. Again, 
the lower borders of the plots correspond to the evaporation 
limits. 

The equilibrium distance r sta b(s) depends on the system's 
parameters as well as on the particle size. Fig. [5] shows solu- 
tions of Eq. ( T29l . computed with and without photophoresis. It 
can be clearly seen that photophoresis significantly increases 
the radii of stable orbits - as expected: photophoresis pushes 
particles away from the star. The curves exhibit a character- 
istic shape: after branching off from the solutions computed 
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without photophoresis, they flatten. While for hLhG there is 
only a small region in which the curve is flatter, the two low- 
luminosity systems have well-developed plateaus, the length 
and distance from the star of which increase with gas den- 
sity. The reason for this is explained below. Note also, that 
the plateaus in the r sta b(s)-curve correspond to areas, where 
solids have a higher surface density. At a certain particle size, 
the curves drop to zero suddenly. This is a consequence of 
mean free path becoming smaller than particle sizes - in this 
area, Eq. (O is not a good approximation anymore, and Eq. ® 
should be used instead to calculate the photophoretic accelera- 
tion. Because in the low Knudsen number regime photophore- 
sis decreases in strength as body size increases, Eq. d28l l and 
( f29T > do not have solutions above a certain body radius, which 
causes r sta b (s) to drop to zero, so that there are no stable orbits 
for solids above this critical size. These considerations are visu- 
alized in Fig. [6] where rj and % are plotted as functions of r for 
different radii s. The third force parameter, j3, can be neglected 
for large particles (s ^> 100/xm), and therefore is not shown. 
We can identify two cases. For small bodies, the j(r)-curve 
has two intersections with rj(r). This means that in a certain 
size range, d28l>-(l29i> have got two solutions. In our work we 
consider only the one farther out from the star, since in transi- 
tional disks we expect most particles to drift inward from larger 
radii, where they are produced by collisions between left-over 
planetesimals in exo-Kuiper belts. Besides, for three out of four 
systems (hGIL is the exception), the inner solution lies in the 
sublimation zone. As s increases further, the ^(r)-curve moves 
downward until it no longer crosses the 7} -curve, and thus no 
solution of d28ll-(|29ll remains. This means that above a certain 
size no stable orbits exist. 
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Fig. 6. The force ratios rj and % as functions of r for different 
particle sizes, in the hLlG system (HR 4796A). The number of 
intersections between the curves corresponds to the number of 
solutions of ( 



Next, we determine how the exact shape of r s ^ D {s) is influ- 
enced by the choice of system parameters. Therefore, in Fig. [7] 
we plot the stability radii for different choices of L*, po and q. 
Note that, when varying in the bottom panel of Fig. [7] we 



change M+ accordingly. Analyzing the plots in Fig. |7j we can 
identify the following dependencies: 

- Higher gas density po moves the curves' middle parts to 
larger r, which is, of course, a consequence of photophore- 
sis being directly dependent on gas pressure. Also, the lim- 
iting particle size beyond which no stable orbits exist de- 
pends almost exclusively on po, while the influence of the 
other parameters is rather weak. 

- The density exponent q, which determines the slope of p (r), 
also controls the slope of r sta b(s): a lower absolute value 
of q produces a steeper curve, as the gas disk spreads to a 
greater distance from the star. 

- Stellar luminosity L* determines the characteristic shape 
of the r sta b(i) curve. While large luminosities produce a 
monotonically decreasing curve, lower (i.e. solar) ones gen- 
erate a plateau at a certain distance from the star, which 
gives rise to a concentration of particles of different sizes 
in that region. 

The mechanism behind the third effect is to be found in the 
denominator of Eq. ©. As we have seen in Fig. [2] in different 
size ranges three different processes reduce the photophoretic 
effect - heat conduction in small particles, radiation for larger 
ones and heat conduction into the surrounding gas for very 
large objects. If the stellar luminosity is low, radiation becomes 
important only for relatively large sizes, while for smaller ob- 
jects only The is relevant in Eq. ©. As Th c is independent of 
size, the resulting curve tends to run parallel to the s-axis, cre- 
ating the plateau. 

In the following, we formulate these considerations quan- 
titatively. We are going to compute solutions of Eq. d29i > 
for very large and very small particles, using the simplified 
large-Knudsen approximation for the photophoretic accelera- 
tion given in Eq. (O. We need to calculate the following ratios: 



Z(Thc) 



and 



-adj 



where a 



flphot,T hl 



■*grav 



flphot.T ra 



phot,T ht 



and a 



phot.T rad 



(30) 



(31) 



denote the accelerations due to 



(large Knudsen regime) photophoretic force using exclusively 
Yhc and Y la d in the denominator, respectively. The first case is 
a good approximation for small particle radii (s < 1 mm, com- 
pare Fig. [2), the second one for larger, but sub-meter objects 
(1 mm < s < lm). 

The quantities xC^hc) and ^(T ra d) can be computed from 
the following expressions: 



Z(Y hc ) = 1.0x10 



-1 



(JhX 
V2.34/ 



Pa 



V lOOOkgirr 



Po 



1AU 



Z(T hc ) 



1AU 



(32) 



and 
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Fig. 7. Dependence of the equilibrium distance r sta t,(.s) on the 
gas density po, the density exponent q, and stellar luminosity 
L*. Top: po running, q and fixed. Middle: q running, po and 
L* fixed. Bottom: L* running, q and po fixed. The parameter 
values are indicated on top of each panel. 
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(33) 



It is easy to show that in the first (small-particle) case, Eq. 
can be solved for s: 



■^stab, small 



B 



(34) 



1/xm E (r[AU])'/2- Z (T hc ) (r[AU])« ' 

For definition of B and E, see Eqs. ( fTTb and ( fTo*b . This allows 
us to calculate the size of particles on stable orbits at a certain 
stellar distance r. Since q is negative, s s tab, small tends to infinity 
as r approaches zero. This produces the curve's plateau, i.e. the 
concentration belt. From Eq. (|34| |. the belt's radius computes 
to: 



/"belt = 



z(T hc ; 



1/(9-1/2) 



(35) 



For the four systems considered, the resulting distances are 
listed in Table [2] 

Table 2. Existence and parameters of the belt. 



System 


rbelt/AU 


C 


1L1G ('older' SN) 


2.0 


0.8 


lLhG ('younger' SN) 


6.5 


2.6 


hLlG (HR 4696A) 


4.0 


0.1 


hLhG (HD 141569A) 


22.8 


0.07 



In the second (large-particle) case, the corresponding ex 
pression 

Sstab.largeM _ B (r [AU])- x ' 2 + Z(T rad ) (r[AU])« +1 



1 fim 



(36) 



can be derived. 

The condition for the formation of the belt is that at its 
position, the value of .s s tab, large has to be sufficiently larger 
than s s tab,smaii- If that is the case, the asymptotic behavior of 
■Sstab, small is seen in the 'exact' curve too, otherwise it is over- 
ridden by J st ab, large ( see Fig- ID- To evaluate this condition nu- 
merically, we set 



r belt — r belt 



Ar 



(37) 



where for Ar a sufficiently small value has to be chosen. We 
use Ar = 0.1 AU. This is necessary because s s tab, small cannot be 
computed at rbelt, as the denominator of Eq. ( f34T > becomes zero 
at this point. Then, the ratio of the two approximating functions 
at rL. lt can be used as a measure of how pronounced the belt is: 



C 



(38) 



•Sstab,large( r belt) / 5 stab, small ( r belt)- 

The values of C are listed in Table|2] Comparing with Fig. [5] we 
find that for C <C 1, the plateau does not appear at all, or is only 
marginal. For C « 1, the belt is well developed. The particle 
size range it encompasses and its degree of concentration (i.e. 
the 'flatness' of the curve in that area) increase with C, which 
makes it a direct measure for the system's tendency to produce 
a belt. In fact, other choices for Ar are possible - for them, the 
values with which C has to be compared, have to be changed. 
For Ar = 0.1 AU, the critical value above which well-formed 
belts appear is C = 1 . 
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Fig. 8. 'Exact' function r sta t, (solid lines) and its two approx- 
imations (dashed lines: Eq. [34l dotted: Eq. l36l l. for the lLhG 
(thick lines) and the hLhG (thin) systems. Note how the be- 
havior of r sta b, small produces the curve's plateau - in the higher 
luminosity case, r sta (, we 'takes over' at smaller particle sizes, 
thus keeping r sta b, small from producing the effect. 



We see that low luminosities and high gas densities are cru- 
cial for belt formation. Since T ra( j °< T 4 « L*, with rising lumi- 
nosity T ra( j becomes important at smaller particle sizes, super- 
seding The before the plateau of s s tab, small is reached. Because 
they increase the overall strength of photophoresis, higher gas 
densities push the plateau outward, reducing its slope - and 
thus increasing the degree of particle concentration in the belt. 

Fig. [8] demonstrates how the exact solution of Eq. (T29b is 
approximated by Eq. d34l > and Eq. d36l > in different particle size 
ranges for the lLhG and hLhG models. 

5.3. Particle size range 

As we have seen, particles of different sizes are influenced 
by photophoresis to a different extent. Its relative importance 
depends on their particle radius, declining for very small 
(micrometer-range) and very large (above meter-range) ob- 
jects. In between, photophoresis plays an important role. In 
order to analyze the critical size range, we plot the ratio x = 
KphotA'r.nophotl as a function of s in Fig. [9] where v r . p hot and 
Vr,nophot are the radial velocity with and without photophoresis 
respectively. The distance is set to r = 10AU. Note that the in- 
terval between the two peaks corresponds to the size range in 
which the two v r have got different signs, because the stability 
radius with photophoresis is greater than 10 AU, while the one 
calculated without photophoresis is smaller. 

If we use the criterion that |jc — 1| > 0.001, photophore- 
sis has to be taken into account in the size ranges listed in 
Table [3] Note, however, that for s < 10/im, the above evalu- 
ation is no longer necessarily valid, as for sma ll particles, the 
photo phoretic force can reverse its direction dTehranian et al. 
l200ll) . Also, for hLhG (HD 141569 A), /3 (s min ) is greater than 
0.5, thus these particles are j3 -meteoroids which escape from 
the system. 
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Fig. 9. Ratio of radial velocities, computed with and without 
photophoresis. The stronger the deviation from unity, the more 
important is photophoresis. The stellar distance is r = 10 AU 

Table 3. Size range of particles affected by photophoresis. 



System 


^min 


Smax 


1L1G ('older' SN) 


9/lm 


3 cm 


ILhG ( 'younger' SN) 


1 /im 


1.5m 


hLIG (HR 4696A) 


lOflm 


3.2cm 


hLhG (HD 141569A) 


1 


6m 



6. Rotation of particles 

Until now, we considered only nonrotating particles. We have 
to check whether this assumption is realistic, since the rotation 
of bodies can transport thermal energy from the dark to the lit 
side, thus banishing the photophoretic effect. There is a number 
of different mechanisms that can induce particle rotation: 

- Collisionally induced rotation. During the formation of 
bodies as well as during their life, they collide with other 
particles, which can change their orbits as well as transfer 
angular momentum to them, spinning them up. 

- Rotation induced by gas drag or radiation forces. If par- 
ticles are not exactly spherical, any force acting on their 
surface will change their angular momentum. For instance, 
radiation pressu re may spin-up the par ticles (the so-called 
windmill effect, Paddac k & Rhedl 19751). but ma y also sta- 
bilize/align then (IDraine & Weingartnerl 1 1996b . Such ef- 
fects are beyond the scope of this paper. 

Whether rotation is able to subdue photophoresis depends on 
four different timescales: f co n, the typical time between two col- 
lisions; f s , the gas coupling time which also determines the time 
needed to slow down rotation; fh ea t> the thermal relaxation time 
needed to establish a stable heat gradient within the particle; 
and f rot , the typical rotation period. 

The collisional time f co u can be estimated in a standard way: 

fcon = — , (39) 

"pVCcoll 
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where n p is the particle number density, v the collisional veloc- 
ity and a co n = 7t(s\ +S2) 2 is the collisional cross section for 
spherical particles of radii s\ and S2- With 



Pd 



Pbuik 5 



(40) 



and the standard assumption Pd/p% = 10 2 , we get for particles 
of equal size (cf. iRrauss et al . 2006, their Eq. 15): 



fcoii — 33 



PbulkS 



(41) 



The stopping time t s was computed in section 13.31 Eq. ( fT9] l. 
Finally, the thermal relaxation time fheat is given by (see 
Krauss & Wurmll2005l Eq. 6): 



fheat : 



PbulkCd^ 



(42) 



where « 1000 Wm^'K" 1 is the particle's heat capacity. 

We can use two different approaches to deal with rota- 
tion. One is to compare collision time to stopping time: if 
f coii A s ^ 1 > damping of rotation occurs faster than excitation by 
impacts, and therefore photophoresis is not significantly weak- 
ened. Using Eqs. (T4TT > and we get 



fcoll 

u 



25 



(43) 



Another option is to check whether rotation, once it has been 
excited by collisions, is sufficiently slow to allow a stable heat 
gradient to be established within the particle: this is the case 
if frotAheat ^> 1. In order to find a rough estimate for t m u we 
assume that the particle transforms the entire kinetic energy of 
its radial motion v r into rotational energy: 



2 L ot ' phere 



1 2 

-mv r 



(44) 



with the solid sphere's moment of inertia /sphere = (2/5)ms 2 , 
we get 



frot 



8 ns 
5 v r 



(45) 



Fig. [10] shows ? C oiiAs as a function of v and f ro tAheat as 
a function of s. The distance is chosen to be r = 10AU. 
We see that collisional timescales are longer than stopping 
timescales by many decades, thus allowing rotation to decline 
between collisions. Nonetheless, rotation frequencies can be 
high enough to cancel out stable temperature gradients. We 
conclude that photophoresis can be suppressed directly after 
collisions through rapid rotation, but remains effective for most 
of the time because gas drag damps spin faster than collisions 
can excite it. 

7. Conclusions and discussion 

In this paper, we have analyzed the effect of photophoresis on 
the dynamics of solid particles in the optically-thin, yet suffi- 
ciently gas-rich, transitional disks around young stars. To this 
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Fig. 10. Ratios of typical timescales. Top: t co n/t s as a function 
of v. Bottom: / rot //heat as a function of s. 



end, we employed a single-body approach and added the pho- 
tophoretic force to the standard array of perturbing forces (stel- 
lar gravity, direct radiation pressure, and gas drag). Particle- 
particle interactions as well as particle growth were not in- 
cluded. We find that photophoresis may lead to noticeable cor- 
rections to the results obtained with models that do not take it 
into account. 

Our main results can be summarized as follows: 

1 . Both with and without photophoresis, solid objects migrate 
inward or outward, until they reach the stability distance 
'stabM, where s is the particle radius. At that distance, ra- 
dial forces cancel each other in the particle's own inertial 
system, and the orbital velocity is equal to that of the gas. 
The stability distance is a decreasing function of s, there- 
fore particles are sorted according to size, with larger bod- 
ies accumulati ng closer to the star. These res ults fully agree 
with those bv lTakeuchi & Artvmowicz](l200ll) . 

2. Photophoresis increases the stability radii, moving objects 
to larger radial distances. The effect is noticeable in the 
size range from several micrometers to several centime- 
ters (for older transitional disks) or even several meters (for 
younger, more gaseous, ones). 

3. The steady-state distribution of solids is completely charac- 
terized by the function r sta b(s), the shape of which depends 
on the system's parameters: 
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- Higher gas densities move the curve to larger distances 
without changing its overall shape. Gas density also 
controls the maximum particle size up to which stable 
orbits can exist. 

- The steepness of the gas density radial profile deter- 
mines the slope of r sta b(s), with flatter profiles gener- 
ating steeper curves. 

- Stellar luminosity determines the curve's shape. While 
high luminosities (~ 2OL ) produce simple decreasing 
curves, low (solar) luminosities generate a plateau at a 
distance from the star that can be computed with the aid 
of Eq. (|35l l. In this area, objects in a certain size range 
may accumulate. 

4. Particle rotation tends to reduce the photophoretic effect. 
Our estimates suggest, however, that it is damped by gas 
drag quickly enough to keep photophoresis at work. 

Our work predicts the formation of a particle concentration 
belt at a certain distance from the star. For the high luminos- 
ity systems, it is not very pronounced (hLhG) or does not ap- 
pear at all (hLlG). Furthermore, the radii of the observed rings 
around HD 141569A are an order of magnitude larger than that 
of the slight concentration belt predicted for the hLhG system. 
It is not likely therefore that the observed structures around 
HR 4796A and HD 141569A are caused by photophoresis. 
While photophoresis is probably active in transitional disks, 
the circumstellar rings of HR 4796A and HD 141569 must 
be shaped by other forces and effects, such as gravitational 
sculpting by planets or interactions with stellar companions. 
Alternatively, a rapid decline of gas density at the disks' edges 
dTakeuchi & Artvmowiczll200ll) or a recently proposed dust- 
gas instability ( Besla & Wull2007l) may cause particles to accu- 
mulate there. 

The model presented here is rather exploratory and rests on 
a number of simplifying assumptions. In the future, we plan to 
investigate the problem a more realistic way, lifting some of 
the assumptions we made to keep the problem tractable. First, 
we plan to deal with particle-particle interactions, taking col- 
lis ions and growth int o account. This can be done in the style 
of iRrauss et al.1 d2006l) . increasing the radius through an expo- 
nential ansatz s = sq exp(f/?o). The latter corresponds to the 
assumption that the object moving through the nebula collects 
smaller particles on its surfa ce. A more detailed appr oach will 
employ statistical methods ( Krivov et al.ll2005l I2006I) . Then, 
we wish to explore rotation in greater detail, calculating colli- 
sion timescales and rotation frequencies using models for two- 
body collisions. Also, the effects of gas drag and radiation 
forces o n non-spherical objects need to be taken into account 
(see e.g. Xu et al. 1999!). Other issues include the variation of 
physical parameters (density, thermal conductivity, emissivity, 
Ji) with size and distance from the star, as well as the global 
evolution and clearing of dust in the system which defines the 
time around which photophoresis can come into play. 

In spite of these unknowns, we have demonstrated that pho- 
tophoretic force in transitional circumstellar disks cannot be 
neglected, and has to be included in elaborate models of such 
systems. 
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